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The antiferromagnetic phase of two-dimensional (2D) and three-dimensional (3-D) Hubbard model 
with nearest neighbors hopping is studied on a bipartite cubic lattice by means of the quantum 
SU(2)xU(l) rotor approach that yields a fully self-consistent treatment of the antiferromagnetic 
state that respects the symmetry properties the model and satisfy the Mermin- Wagner theorem. 
The collective variables for charge and spin are isolated in the form of the space-time fluctuating 
U(l) phase field and rotating spin quantization axis governed by the SU(2) symmetry, respectively. 
As a result interacting electrons appear as a composite objects consisting of bare fermions with 
attached U(l) and SU(2) gauge fields. An effective action consisting of a spin-charge rotor and a 
fermionic fields is derived as a function of the Coulomb repulsion U and hopping parameter t. At 
zero temperature, our theory describes the evolution from a Slater {U <g t) to a Mott-Heisenberg 
{U S> t) antiferromagnet. The results for zero-temperature sublatice magnetization (2D) and finite 
temperature (3D) phase diagram of the antiferromagnetic Hubbard model as a function of the 
crossover parameter U /t are presented and the role of the spin Berry phase in the interaction driven 
crossover is analyzed. 



I. INTRODUCTION 

Effective low-energy theories are frequently employed 
in essentially all fields of physics. In the realm of strongly 
correlated electrons, spin only Hamiltonians, are exam- 
ples of effective low-energy theories that apply in the 
limit of strong interactions. Due to the numerical ex- 
pense in solving these models for large lattices, it is im- 
perative to look for methods which relay on advanced an- 
alytical approaches. A key question in this context con- 
cerns the emergence of low energy scales, much smaller 
than the bare Coulomb interactions between the elec- 
trons, which govern the existence and the competition of 
different phases. This can be studied by considering pro- 
totypical lattice models of strongly correlated electrons. 
In the strong coupling limit, the t — J Hamiltonianfi de- 
rived from the large— f7 Hubbard model;^ is often used to 
describe the low lying excitations. At half-filling, it di- 
rectly reduces to the quantum antiferromagnetic Heisen- 
berg model. However, for finite U as in the the Hub- 
bard model, the quantum-mechanical objects are not lo- 
cal spins since one has still mobile electrons and one ex- 
pects that calculation of the ground-state phase diagram 
as a function of the density and interaction strength is 
more difficult than for the Heisenberg model. In the weak 
coupling limit (?7 <C i), a Fermi surface instability gives 
rise to a spin-density-wave ground state as described by 
Slater)^ where the antiferromagnetic (AF) long-range or- 
der produces a gap in the quasiparticle excitation spec- 
trum. In the strong-coupling regime (U ^ t), fermions 
are effected by the strong Coulomb repulsion causing 
the Mott-Hubbard localization.'* The created local mag- 
netic moments subsequently order at low temperature 
and give rise to a Mott-Heisenberg antiferromagnet. It is 
well known that even for small but finite interactions the 



Hartree-Fock (HF) Neel temperature is proportional to 
U, which is unrealistically high since C/ is a large energy 
scale, of the order of electron volts. This wrong predic- 
tion had to be expected since correlations are absent in 
the HF approach. The requirement of self-consistency by 
incorporating spin and charge fluctuations, while main- 
taining the essential spin-rotation symmetry, summarizes 
the challenging nature of magnetic ordering in strongly 
correlated systems. A variety of theoretical methods are 
available for the study of strongly correlated systems. In 
weak or strong coupling perturbative treatments usually 
are usedi^ On the other hand numerical - very effective in 
low dimensional models - are limited to finite systems, 
which requires an extrapolation to the thermodynamic 
limit that is often problematic. As a result the numerical 
approach docs not provide in general a unifying picture, 
which only analytical approaches can give. 

The objective of the present paper is to quanti- 
tatively investigate correlation effects in a antiferro- 
magnetic state of the Hubbard model within a spin- 
rotationally-symmetric scheme that is fully compatible 
with the Mermin- Wagner theoremJi To this end, we de- 
scribe a theoretical approach which provides a unified 
view of two {2D) and three (3D) dimensional model at 
half-filled Hubbard model for any value of the Coulomb 
repulsion U, which is able to handle the evolution from 
the Slater to the Mott-Heisenberg antiferromagnet that 
captures correctly both the spin and charge degrees of 
freedom. We address the above questions by implement- 
ing the charge-U(l) and spin-SU(2) rotationally invari- 
ant handling of the Hubbard model. By recognizing spin 
and charge symmetries we explicitly factorizc the charge 
and spin contribution to the original electron operator 
in terms of the corresponding gauge fields that leads to 
a composite particle, which is the union of an electron 
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with U(l) and SU(2) gauge potentials. In this scheme 
the charge and spin excitations emerge in terms of a 
U(l) phase and variable spin quantization axis: the ef- 
fective field theory for the strongly correlated problem is 
thus chracterized by the U(2)=U(l)xSU(2) group, where 
the gauge potential in U(l) describes the evolution of a 
particle scalar characteristic, which is naturally associ- 
ated with an electric charge, while the gauge potential in 
SU(2) describes the nontrivial dynamics associated with 
the evolution of the vector internal characteristic of a 
particle such as spin. 

The outline of the paper is as follows: In Sees. II we 
introduce the model and, in Sec.III-V, we develop the an- 
alytical background needed for the calculations. In Sec- 
tion VI we find a closed set of self-consistent equations 
for the antiferromagnetic gap and order parameter, while 
the numerical evaluation of self-consistent equations is 
presented in Sees. VII, where the phase diagrams con- 
sidering the paramagnetic the antiferromagnetic phase 
for the Hubbard model and in different dimensionality 
are calculated. 



II. THE MODEL 

Our starting point is the purely fermionic Hubbard 
Hamiltonian Ti. = Tit + TLu'- 



n 



-t 



E 



[4(r)c„(r') +h.c.] + C/^ n|(r)nx(r) 



(1) 



Here, (r,r') runs over the nearest-neighbor (n.n.) sites, 
t is the hopping amplitude, U stands for the Coulomb 
repulsion, while the operator cj,(r) creates an electron 
with spin a =^,i at the lattice site r, where ric((r) = 
cjj(r)cct(r). Usually, working in the grand canonical en- 
semble a term is added to Ti. in Eq. ([T]) to control the 
average number of electrons, H ^ H — fJ, J2r '^('^) '^'^^ 
fi being the chemical potential and 7i(r) = n^{r) + rix(r) 
the fermionic number operator. 



A. Grassman action 

The functional integral representation of models for 
correlated electrons allows us to implement efficiently 
the method of treatment. It is customary to introduce 
Grassmann fields, CQ,(rr) depending on the "imaginary 
time" < T < P = l/ksT, (with T being the temper- 
ature) that satisfy the anti-periodic condition Ca{rT) — 
—Ca{rT + (3), to write the path integral for the statistical 
sum Z = J [VcDc] e~'^[^''^l with the fermionic action 



S[c,c]=Sb[c,c]+ / dTT-l[c,c] 
Jo 

that contains the fermionic Berry termii 

Sb[c,c] dTCa{rT)drCairT) 
v„ "'0 



(2) 



(3) 



which will play an important role in our considerations. 



III. SU(2)xU(l) ACTION 

It is customary to introduce auxiliary fields for the spin 
and charge fluctuations via a Hubbard- Stratonovitch 
(HS) transformation to decouple the interaction term in 
the Hubbard Hamiltonian. However, such a procedure 
usually leads to a loss of the spin rotational invariance. 
For strongly correlated system in order to properly ac- 
count for the nature of elementary excitation it is crucial 
to construct a formulation of the theory which naturally 
preserves the existing symmetry present in the Hubbard 
Hamiltonian. For this purpose the density-density prod- 
uct in Eq. ([T]) we write, following Ref. iSs, in a spin- 
rotational invariant way: 



(4) 



where S"'(rT) = ^J2aa'^ai'^^)^aa'^c'{rT) dcuotcs the 
vector spin operator (a = x, y, z) with ct" being the Pauli 
matrices. The unit vector 

f2(rr) = [sin'!9(rT) cos c/?(rr), sini9(rT) sin iy9(rT), 

cos?9(rT)] (5) 

written in terms of polar angles labels varying in space- 
time spin quantization axis. In order to maintain 
spin rotational invariance, one should consider the spin- 
quantization axis to be a priori arbitrary and integrate 
over all possible directions of r2(rT) in the partition fimc- 
tion. The reason for choosing the decomposition in Eq. 
(131) is that it allows to recover the Hartree-Fock solution 
at the saddle point level within the functional integral 
formalism. By decoupling spin and charge density terms 
in Eq. ([4]) using auxiliary fields g{rT) and iV{rT), re- 
spectively we write down the partition function in the 
form 

Z = y" [Vn] J [VVVq] j [VcDc] X 



X e 



-S[fn,V,g,c,c] 



(6) 



where [2?n] = Ilrr, is the spin- 

angular integration measure. The effective action reads: 



5[JI,l/,e,c,c] = E / 



U 



U 



iV{rT)n{vT) + 2£»(rr)f2(rr) • S(rT)] 

(7) 



-I- Sb [c, c]+ [ drTLt [c, c] . 
Jo 



We devise a systematic way of decomposing the fluctuat- 
ing fields contained in the action in Eq. ([7]) that enables 
us to obtain a low energy effective theory. In the follow- 
ing introduce bosonic fields describing charge and spin 
fluctuations, and the fluctuating spin-quantization axis. 
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A. U(l) rotor charge frame 



B. Rotating SU(2) spin reference frame 



We observe now that the Hubbard Hamihonian has 
a local U(l) gauge symmetry, when expressed in terms 
of the underlying electron variables. This points out a 
possibility of an emergent U(l) dynamical gauge field 
as a fluctuating complex field attached to fermionic 
variables, which is dynamically generated, by interact- 
ing fermions. In the modern language it is called a 
fermion-flux composite. Technically, the appearance of 
the U(l) field is based on the Hubbard-Stratonovich de- 
coupling of the four-fermion interaction - a typical way 
to "bosonize" a fermionic system in higher-dimensional 
problems. The essence of the method is to eliminate a 
mixed fermion-boson term in the resulting action by a 
gauge transformation!^ The U(l) formulation begins by 
rewriting the electron as the product of a charge-neutral 
fermionic spinon and a charge-carrying U(l) rotor, which 
is constrained to lie on the unit circle in the complex 
plane. To this end, we write the fluctuating "imaginary 
chemical potential" iV{rT) as a sum of a static Vb(r) and 
periodic function V{rT) — Vb(r) + V^(rr) using Fourier 
series 



-. oo 



(8) 



with ujn = 2im/l3 {n = 0, ±1,±2) being the (Bose) Mat- 
subara frequencies. Now, we introduce the U(l) phase 
field (/>(rr) via the Faraday-type relation 

OT i OT 

Furthermore, by performing the local gauge transforma- 
tion to the new fermionic variables fa{rT): 



Ca 


(rr)- 




Ca 


(rr)_ 





z(rr) 
z(rT) 



Urr) 
Uvt) 



(10) 



where the unimodular parameter |z(rr)p = 1 satis- 
fies z(rr) — e*'^'-""^'', we remove the imaginary term 
i drV {vT)n{vT) for all the Fourier modes of the V{yt) 
field, except for the zero frequency. Accordingly, the in- 
tegration measure over the group manifold becomes 



{m(r)} 



0(r/3)=0o(r)+2irm(r) 



d0o(r) / Vcj^iyr). 



(rO)=0o(r) 



(11) 

Since the homotopy group 7ri[[/(l)] forms a set of inte- 
gers, discrete configurations of (/"(rr) matter, for which 
(j){r(3) — 4>{rO) = 27rm(r), where rn(r) = 0, ±1,±2,.... 
Here, to e Z labels equivalence classes of homotopically 
connected paths. Thus the paths can be divided into 
topologically distinct classes, characterized by a winding 
number defined as the net number of times the world 
line wraps around the system in the "imaginary time" 
direction.— 



Subsequent SU(2) transformation from /a(rr) to 
ha{vT) operators. 



/i(rr) 
/2(rr) 



R(rr) 



h2(vT) 



(12) 



takes away the rotational dependence on Jl(rT) in the 
spin sector. This parametrization makes clear that the 
SU(2) matrix rotor is identical to the more familiar 0(4) 
rotor, a quantum particle constrained to the 3-sphere 



R(rr) 



cos (I) 
sin (I) 



eiiv+x) cos 



sin (I) 
I) 



(13) 



with the Euler angular variables V3(rr), i?(rr) and x(rT), 
respectively. The link between S0(3)and SU(2) rotations 
is established by means of the Hopf mapii 



R(rr)6-^Rt(rr) = & ■ n{rT) 



(14) 



that is based on the enlargement from two-sphere S2 to 
the three-sphere ^3 ^ SU{2). Here, the extra variable 
x(rT) represents the U(l) gauge freedom of the theory as 
a consequence of S2 — > S3 mapping. One can summarize 
Eqs and (fT^ by the single joint gauge transforma- 
tion exhibiting electron operator factorization 



where 



U{vt) = z(rr)R(rT) 



(15) 



(16) 



is a U(2) matrix which rotates the spin-charge quanti- 
zation axis at site r and time t. Eq. (fT5|) reflects the 
composite nature of the interacting electron formed from 
bosonic spinorial and charge degrees of freedom given 
by i?QC['(rT) and z(rT), respectively as well as remaining 
fermionic part ha (rT) . In the new variables the action in 
Eq. ([7]) assumes the form 

S\p,,4',Q,h,h\ ~ SB[h,h]+ I dTT-ln,4,[p,h,h] 

Jq 



-^5o[<^]+2V / dTA(rr) • S^(rT), 



(17) 



where S;j(rr) = ha{TT)&ajhj{rT). Furthermore, 



„ "'0 



0'(rT) 
U 



1 2fi 
111' 



(18) 



stands for the kinetic and Berry term of the U(l) phase 
field in the charge sector. The SU(2) gauge transforma- 
tion in Eq. and the fermionic Berry term in Eq. ^ 
generate SU(2) potentials given by 



d 



d 



d 



K\vT)drK{YT) = RT ( + ^_ + ) R 



-a-- A(rT), 



(19) 
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where 



IV. SPIN-ANGULAR ACTION 





2 ^ ' 


sinx(rT) 




- 




sa\9{TT) cosx(rT) 




Ay[rT) = 


2 ^ ' 


cosx(rT) 




+ 




sin 9{yt) sin x(rT) 








cosi?(rr) + -x(rT) 


(20) 



are the exphcit expression for the vector potential in 
terms of the Euler angles. 



C. Integration over Vb and g 



Since we are interested in the magnetic properties of 
the system a natural step is to obtain the effective ac- 
tion that involves the spin-directional degrees of freedom 
f2, which important fluctuations correspond to rotations. 
This can be done by integrating out fermions: 



[V(t>Vn] I [DhVh] e 



-S[ip,<ti,-d,h,h] 



where 



= - In / [V(j)VhVh 



(24) 



(25) 



generates the cumulant expansions for the low energy 
action in the form S[Vl\ ~ Sb[^] + 5,/[r2]. 



We observe that the spatial and temporal fluctuations 
of the fields Vb(r) and q{yt) will be energetically penal- 
ized, since they are gapped and decouple from the an- 
gular and phase variables. Therefore, in order to make 
further progress we subject the functional to a saddle 
point HF analysis: the expectation value of the static 
(zero frequency) part of the fluctuating electrochemical 
potential Vo(?') we calculate by the saddle point method 
to give 



K,(r) 



U 



(21) 



where n = J2a(^a{YT)ha{rT)) and the saddle point with 
respect to p gives 



p(rr) - (-l)'-A, 



(22) 



with Ac setting the magnitude for the Mott-charge gap. 
The choice delineated in Eq. corresponds to the 

saddle point of the antiferomagnetic (with staggering Ac) 
type. Note, that the notion "antiferomagnetic" here does 
not mean an actual long-range ordering - for this the 
angular spin-quantization variables have to be ordered 
as well. The fermionic sector, in turn, is governed by the 
effective Hamiltonian 

r 

{r,r'> 

~p,'Yha{rT)ha{rT), (23) 



where fl = fi — nU /2 is the chemical potential with a 
Hartree shift originating from the saddle-point value of 
the static variable Vq. 



A. Topological theta term 

In general, in addition to the usual exchange term, 
the action describing antiferromagnetic spin systems is 
expected to have a topological Berry phase term 



Sein] = -2E / dTAivT) ■ {Sh 
, Jo 



(r'r')>, (26) 



(5^(rr)) = i(nT-nx) = ^, (27) 



where 



which results from the saddle point value in Eq. (|22p . In 
terms of angular variables, the Berry term becomes 

5b = - V / dT [^(rr) cos ^(rr) + xirr)] . (28) 

If we work in Dirac "north pole" gauge xi^''') = 
— i^(rr) one recovers the familiar form 5s [fi] — 
1 Sr lo dTip{YT) [1 — cos ?9(rT)] . Here, the integral of the 
first term in Eq. ([28]) has a simple geometrical inter- 
pretation as it is equal to a solid angle swept by a unit 
vector r2('i?, Lp) during its motion. The extra phase factor 
coming from the Berry phase, requires some little extra 
care, since it will induce quantum mechanical phase in- 
terference between configurations. In regard to the non- 
perturbative effects, we realized the presence of an addi- 
tional parameter with the topological angle or so-called 
theta term 

Ac 

6 = ^ (29) 

that is related to the Mott gap. In the large-?/ limit 
one has Ac so that Q ^ \ relevant for the half- 

integer spin. However, for arbitrary U the theta term 
will be different from that value, which, as we show will 
be instrumental for destruction of the antiferomagnetic 
order away from the spin- localized U ^ oo limit. 
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B. AF exchange second order cumulant term in the hopping element t 

containing four fermion operators: 

Now we proceed with the calculation of the exchange 
term in the spin-angular action. We concentrate on the 



|ri — r, I — n.n 



X ^(r2r')z(r'2r')^[Rt(r2r')R(r'2r')]^^,/i^(r2r')V(rir')), 

|r2— r2l=n.n. 77' / 



(30) 



where (...) denotes averaging over U(l) phase fields and 
and fermions. The averaging in the charge sector is per- 
formed with the use of the U(l) phase action in Eq. ^ 
to give 



X exp 



u 



[r-r'f 

/3 



(31) 



Furthermore, with the help of the gradient expansion 

R(rr') = R(rr) + (r' - T)9^R(rT) 

+ 0[(r'-r)2] 

h(rT') = h(rT) + (t - T)arh(rT) 

+ 0[{T'-Tf] (32) 

we write the relevant part of the action in the form 

r^^ \ U 1 

Sj[VL] = -- drdr'exp --|r-T'| 
^ Jo L ^ 

X E E[l^^(r^)R(^'^)]o«'[R^(r'r)R(rr)]^^, 



|r— r' |— a 

-fi' 

x(ha{vT)hY{vT)){ha'{v'T)h^{v'T)). 



(33) 



In the low temperature limit (on the energy scale given 
by [/), by making use of the formula 



/3 

lim / dr' e 



\t-t'\U 



au 

2 2e-T- 



(34) 



/o U U 

and with the aid of the fermionic occupation numbers 
(/ia'(r'r)/i^(r'T)) = (1 - 7i^)5a''i 

{ha{rT)hy{rT')) ^ UaSaY (35) 

we arrive at 

ir rRtfrrlRfr'rV, 



= \j j dr [RUrr)R(r'r) 

|r-r'|=a "7 

X [Rt(r'r)R(rr)] (n^ - l)n„. (36) 



Finally, making use of the following composition formula 
for the SU(2) matrices 



[Rt(rr)R(r'r)]^^ [Rt(r'r)R(rr)] 
1 



a7 L V ' ^ ' J 70: 

[l-f2(rr).fl(r'r)](l-<5„^) 



+ i[l + f^(rT) •0(r'T)],5„^ 
we obtain the desired part of the spin action 

(rr') " 

-I- n(n-2)] 
with the AF-exchange coefficient 



(37) 



(38) 



J(A, 



The factor ~ comes from integration over U{1) 

charge degrees of freedom (see, Eqs. (pij) . ([M)) ). whereas 
the occupation numbers result from integration over 
fermionic variables (see, Eq. Wf\ . From the Eq. ([551) 
it is evident that for f/ ^ 00 one has J(Ac) ~ ^ 
since ^ 1 in this limit. In general the AF-exchange 
parameter persists as long as the charge gap Ac exists. 
However, J(Ac) diminishes rapidly in the U /t — s- weak 
couphng limit, see Fig. [TJ 



V. FERMIONIC SECTOR 

Now we evaluate the effective interaction between 
fermions by integrating out by means of cumulant ex- 
pansion the gauge degrees of freedom. To this end we 
write the partition function as 



[V(j}DVl\ / [DhVh] 



-S[ip,it>,S,h,h] 



[VhVh] 



'S[h,h] 



(40) 



6 




FIG. 1: (Color online) The antiferomagnetic exchange param- 
eter J as a function of the Coulomb interaction U in two and 
three dimensions. The fermionic occupation number is fixed 
by n = 1. 



where 



S[h, h] = - In / [D(f>Vn]e 



^-Slifi,<t>,^,h,h] 



(41) 



generates cumulant expansion for the effective fermionic 
action. Since, both U(l) and SU(2) gauge fields couple 
to the hopping element t, in the lowest order of cumulant 
expansion we reveal the hopping renormalization t tg, 

9 = 9c9s = 

5c = {z{yt)z{v't)) 

5. = ([Rt(rr)R(r'r)]^P 

= ([Rt(rr)R(r'r)]u) (42) 

where the coefficients gc and contribute to the band 
renormalization in a way, which is similar to the action 
of the Gutzwiller factorsji^ and have to be calculated 
self consistently, according to the Eq. that involve 
charge and spin-angular correlation functions. However, 
as long as there is no ordering in the charge sector, 
{z{yt)z{y' t)) = and g = resulting in renormalized 
hopping t = 0. Thus, we have to rest on the second or- 
der of the cumulant expansion, in which one obtains a 
contribution to the effective action in the form 

S^^\h,h]=-— dr V.Ft(rrr'T).F(rrr'T) (43) 
Jo I 

(rr'> 

with the bond operators 

^(rrr'r) = H^^)^^ ^) + H^^)^ ii^' ^) (44) 

^/2 




FIG. 2: (Color online) The effective hopping parameter tj/t 
(see, Eq. (|49ll) related to the AF-exchange as a function of 
t/U at n = 1 for two and three dimensional Hubbard model 
at n = 1. 



Since the iS*-^-* \h, h] is quartic in the fermionic variables, 
we resort to the Hubbard-Stratonovich decoupling with 
the help of the complex variables defined on the links of 
the lattice 



= / [V'^^]e <--'> » , (45) 



where ^ n<rr')r d'^S.irrr'r) and = dRe^dlm^. 
Saddle point with respect to ^ gives 

e= ^(.F(rrr'r)) = ^(/,4rr)/.4r'r)). (46) 



Denoting 



V = y^(/ta(rr)/ta(rV)), 



(47) 



which plays the role of the kinetic energy parameter for 
fermions we obtain 



S^^'^[h,h]^^^+tj ^ [/i„(rr)/i„(r'T) + /i.c.], (48) 

{rr')a 



where 



t. 



-Jv 



(49) 



is the effective hopping parameter that involves the an- 
tiferromagnetic exchange parameter J. The interaction 
dependence of this quantity is depicted in Fig. [2j Sum- 
marizing the results we obtain for the fermionic action 



S[h,h]^SB[h,h]+ dTH[h,hl 
Jo 



(50) 
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with the effective Hamiltonian 



r 

(r,r' ,a) 



(51) 



The result of the gauge transformations is that we have 
managed to cast the strongly correlated problem into a 
system of mutually non-interacting fermions, submerged 
in the bath of strongly fluctuating U(l) and SU(2) fields, 
whose dynamics is governed by the energy scale set by 
the Coulomb interaction U coupled to fermions via hop- 
ping term and with the Zeeman-type contribution with 
the massive field g{rT) related to the Mott gap A^. To 
calculate the latter one has to introduce two inequivalent 
sublattices, let say A and B and write the Hamiltonian in 
terms of two sublattice operators in the reduced Brillouin 
zone (RBZ). The fermionic propagator then reads 



G(ku;„) = 



it^n + ek - A* 
{ivn - ek - p.){iiy„ -f Ck + /l) - 



(52) 



where i/„ = 7r(2n-|-l)//3, n = ±1, ±2, . . . are the fermionic 
Matsubara frequencies. The self-consistency equations 
for the Mott gap, the kinetic energy bond parameter and 
fermionic occupation number are given by 

^= = ^ E 'GTi(k^«) 



— y 

^^kt^. 



1 ' 



(53) 



where the sums with prime index denote the summa- 
tions over wave vectors inside the the RBZ and 7(k) — 
(cos fca; -|- cosky)/2. By performing the summations over 
Matsubara frequencies one obtains explicitly with the use 
of the fermionic distribution np{x) = l/(e'^^ + 1)- 



27V ^ 



(54) 



for the gap parameter, 

2N ^ ' El, ' V y 

k 

the fermion kinetic energy parameter, 

n ^ ^y^inp {-El, - p.)+ rip (E'k - /!)] (56) 



and the occupation number, respectively. At half-filling, 
solutions of the Eq. ((5i|) for the gap are stabilized 
for any arbitrarily small U, however by computing of the 
free energy one can find that no stable antiferromagnetic 
solutions away from n — 1 exists. ^'^ Solutions of the self- 
consistency equation for the gap Ac leading to densities 
away from half-filling correspond to maxima instead of 
minima in the free energy. If a certain occupation near 
the half-filling is enforced, the system will experience a 
phase separation in a half-filled antiferromagnetic and a 
non-half-filled paramagnetic region. 



VI. FORMULATION 

Since the fermionic field can be systematically inte- 
grated out, the main practical difficulty comes from the 
dynamics of spin-directional fluctuations. To proceed 
with the spin-bosonic action we resort to the CP^ repre- 
sentation (see, e.g. Refill. In the CP"'^ representation, 
the SU(2) rotation matrix is expressed in terms of two 
Schwinger bosons. 



R(rr) 



Ci(rr) -C2(rr) 
C2(rr) Ci(rr) 



(57) 



with the constraint |(i(rT)p + |^2(r''')P = 1- The uni- 
modular constraint can be resolved by using the Euler 
angles parametrization 



Ci(rT) 
C2(rT) 



g-fb(rr)-|-x(rr)]^Qg 



1?(rT) 



1?(rT) 



(58) 



which make link between the Ci(rT), C2(rT) fields and 
r2(rT) variables. By definition 



(59) 



a7 



are the "bosonic" spins in the complex-pro jective (CP""^) 
formulation, while the action, see Eq. ([55]) becomes 



I i\ Jo 

(rr'> 



Sc(rr).Sc(r'r)-i 



(60) 



Consequently,_the complete spin-bosonic action C] = 
SeiCC] +Sj[CX] reads 

SeiCX] = -2ey2{-lY / dTUrr)Urr) 
SjICC] = ~jy2 dTA{rTr'T)A{rTr'T), (61) 

{rr') 

where 9 = ^ is the "theta angle" parameter in the 
Berry-phase term, while the AF-exchange term Sj [(, (] 
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wc write with the help of the valence-bond operators 
^(rrr'r) for which the following relations hold 

Sc(rr) • Sc(rV) = -^(rTr'r)>l(rrr'r) + ^ 
^(,,,V)^CTMa(rV)-_a(rr)CT(rV)_ ^^^^ 

A. HS decoupling 

In order to achieve a consistent representation of the 
underlying antiferromagnetic structure, it is unavoidable 
to explicitly split the degrees of freedom according to 
their location on sublattice A or B. Since the lattice is 
bipartite allowing one to make the unitary transforma- 
tion 



Ci(rT) ^ CT(i-r) 
for sites on one sublattice, so that 



(63) 



Furthermore, one then performs a saddle-point approxi- 
mation over the Q and A fields 



ap(rrr'r) = ^ (C„(rr)Ca(r'r)) 

a=l 



= - — (^'(rrr'r)) 
\ = jE(Ca(rr)Ca(rr)) 



(69) 



by assuming the uniform solution (5sp(rTr'T) = Q we 
obtain for the Hamiltonian in the spin-bosonic sector 



Wq[C,C] = ^5]^[Ca(k,a;„),Ca(-k,-a;„)] 

\<.uJn a— 1 



Cc«(-k, -w„) 



(70) 



with 



A{vtv't) ^ A'irTr'T) = J2^^^^^^^p^- (64) 

a— 1 ^ 

Biquadratic (four-variable) terms in the Lagrangian can- 
not be readily integrated in the path integral. Introduc- 
ing a complex variable for each bond that depends on 
"imaginary time" Q(rTr'T) we decouple the four-variable 
terms ^'(rrr'r)^'(rrr'T) using the formula 



J2 J dT{^\Q\''+QC<+QC<) 



<"') 



(65) 



where V^Q = U{rr')r d^Qi^^^''^) and d^Q = 
dReQdlmQ. To handle the unimodularity condition one 
introduces Lagrange multipliers A(rT) at each time and 
site. Then with the help of the Dirac-delta functional 



2i9u!n + A -27k(3 
-^7kQ -2i6u>„ + A 



(71) 



Subsequently, performing the sums over Matsubara fre- 
quencies one obtains 



N ^ 26 2a;k 

k 



1 I ^ \ \ 

k 



46* 



(72) 



where Wk = -\/A^ — (-27k<5)^ and z is the lattice coordi- 
nation number. 



VII. AF LONG-RANGE ORDER PARAMETER 



5(|C(rr)|2-l) = 







/ 


2m _ 



J2J drAdCI^-l) 

e ■■ (66) 



the variables Ci(rT),C2(rr) are now unconstrained 
bosonic fields. Thus, the local constraints are reintro- 
duced into the theory through the dynamical fluctuations 
of the auxiliary A field 



Z = 
where 



[C,C]) 



(67) 



Wq[C,C] = E/ dT{QC-C + QC-C + K-C)- (68) 



A characteristic property of strongly correlated sys- 
tems is the existence of local moments. Weak-coupling 
theories usually fail to properly describe these local 
moments. When electron correlation effects become 
stronger, in general, spin fluctuations have to be con- 
sidered seriously. The HF transition temperature bears 
a physical meaning as a temperature below which the 
amplitude Ac of the AF order parameter takes a well- 
defined value. This is also interpreted as the appearance 
of local moments. However, a nonzero value of A,; does 
not imply the existence of AF long-range order. For this 
the angular degrees of freedom n(rT) have also to be or- 
dered, whose low-lying excitations are in the form of spin 
waves. In the CP^ representation (where the Neel field 
is represented by two Schwinger bosons) Bose-Einstein 
condensation of the Schwinger bosons at zero tempera- 
ture signals the appearance of AF long-range order. The 
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AF order parameter in terms of the original fermion op- 
erators is defined as 

lUAF = ^(-l)'-(5^(rr)) = 

r 

= ^(-fr(n(rr)).(S,(rr)). (73) 

r 

Owing the fact that (5JJ(rT)) = (— l)'"Ac5a,z we obtain 
ruAF = Ac^(l]^(rr)) 

r 

= A,^ [(CT(rr)CT(rr)) - (Ci(rr)Ci(rr)>] (.74) 

r 

Furthermore, the order parameter for the CP^ "boson 
condensate" is 

(CQ(kw„)) = (Ca(ktJ„)) 



I3N 



'Tlo<5o,w„^T,a (^k,0 + (^Icq) • (75) 



This yields a macroscopic contribution (i.e., order one) 
to the staggered magnetization and represents a macro- 
scopic contribution to the CP^ bosons density, of the 
a =] bosons at the mode with k = 0, cij„ = thus giving 



A, 



mAF 



(3UN 



[(CT(kc^«)CT(k^n)) 



-(Ci(kc^„)Ci(kc.„))] = -^ml (76) 

Finally, the fraction of condensed Schwinger bosons is 
given by 

which represents the extension of the saddle point equa- 
tion for the Lagrange multiplier to the region of the or- 
dered state. 



A. d = 2 Hubbard model 

In two dimensions, we expect no long-range AF or- 
der at finite temperatures due to Mermin and Wagner's 
theorem. This could be verified by explicitly performing 
two-dimensional momentum summations in Eqs. (j72p 
with the help of density of states for the two-dimensional 
square lattice P2d{f) = /_^^['^^k/(27r)^(5[e — e(k)], where 
e(k) = Co (cos kx + cos ky) with 



P2.(e) = ( 1 - 



X K 



2^ 



(78) 




2.0 



UA 



FIG. 3: (Color online) The Mott gap Ac, fraction of con- 
densed Schwinger bosons mo and AF order parameter mAF 
for the half-filled Hubbard two-dimensional model at zero 
temperature. 



As a result, there is no antiferromagnetic order at finite 
temperature in two-dimensional Hubbard model. Tak- 
ing the zero temperature limit in Eq. (j72p and fix- 
ing the fermionic occupation number at n = 1, in the 
limit U/t — > (X) we find the order parameter value, 
mAF = 0.308 in the ground state (in agreement with 
the calculations from the Ref. [isl ) that is less than the 
classical value S = 1/2. Monte Carlo calculation on 2D 
Hubbard modeli^ gave mAF = 0.4 This effect is due 
to the quantum zero-point motion, which has its ori- 
gin in the noncommutability of the Hamiltonian and the 
staggered magnetization. In the opposite weak coupling 
limit, U/t — *■ 0, the gap Ac persists at arbitrary small 
value of U/t^ however the true order parameter mAF^ 
which involves also the density of condensed Schwinger 
bosons vanishes sXU/t k. 0.621 (see. Fig. The de- 
struction of the AF order is due to the Berry phase term 
whose coefficient, cf. Eq. ((29|) differs from the local- 
ized spin value S" = 1/2 while entering the weak coupling 
limit. In particular, for U/t 0, that is in the weak cou- 
pling limit, 9 goes to zero along with the charge gap Ac 
(see, Eq. ((29)) ) and the self consistency Eqs. ([72)1 predict 
vanishing of the long-range AF order. 



B. d = 3 Hubbard model 

In three dimensions, for a system with an ordered 
ground state, thermally excited states reduce the spin 
correlations at finite temperatures. When the tem- 
perature is much higher than the typical coupling en- 
ergy scale J, we expect the spins to be uncorrelated 
at large distances and the magnetization mAF to van- 
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0.5 1.0 1.5 2.0 



UA 

FIG. 4; (Color online) Same as in Fig. [2] but for three- 
dimensional Hubbard model. 

ish in the absence of an ordering field. This requires 
a phase transition at some temperature between the 
ordered and disordered phases. As in the previous 
case we employ the density of states for the cubic lat- 
tice p3d(e) = Z^m{2v:)H[e - 6(k)], where e(k) = 
eo (cos kx + cos ky + cos kz ) : 




with oi = max(— 1, —2 — e/eo), 12 — min(l, 2 — e/eo). 

The interaction dependence of the AF magnetic mo- 
ment is depicted in Fig. |31 In the U ^ 00 localized 
limit it is mAF ~ 0.422, i.e less the mean field value 
l^c/U = 1/2, however bigger than in the case of the 
2D Hubbard model. Finally, Fig. [5] displays the cal- 
culated antiferromagnetic phase diagram as a function 
of temperature and interaction strength. At weak cou- 
pling our theory clearly describes a Slater antiferromag- 
net with an exponentially small AF gap. As U increases, 
the Slater antiferromagnet progressively evolves into a 
Mott-Heisenberg antiferromagnet with an AF gap of or- 
der U . In the weak interaction limit there is a destruc- 
tion of the AF order at U/t = 0.676 (see Fig. |6]), due 
to the topological Berry phase term whose coefhcient, cf. 
Eq. deviates from the from the localized spin value 
5 = 1/2 in the weak coupling limit U/t. The AF critical 
temperature has a maximum at U/t « 3.78. It is worth- 
while to compare our results with the previous work on 
the subject. Numerical methods such as dynamical clus- 
ter approximation (DCA)i'I'i-'"' give U/t k, 7.5, whereas 




5 10 15 20 



U/t 

FIG. 5: (Color online) The temperature-interaction phase di- 
agram for the three dimensional Hubbard model at half filling. 
Depicted is the temperature Tg for the vanishing of the gap 
parameter Ac as well as the true critical temperature Tc at 
which the log-range AF order ceases to exist, signalled by 
vanishing of mAF- 

dynamical mean-field theory approximation (DMFA)^^ 
predicts U/t ~ 10. The methods based on a perturba- 
tion theory with respect to the interaction strengt h^°i^^ 
are unable to reproduce the maximum in the AF criti- 
cal temperature as a function of U/t. The significantly 
higher values of U / 1 resulting from DCA and DFMA have 
to be explained by the restricted ability of these methods 
while handling spatial fluctuations. Regarding the value 
of maximum of the critical temperature Tc/t w 0.667 
found here, it agrees with the result of Monte Carlo sim- 
ulations by Scalettar et aP^ Tc/t « 0.72 and by Hirsch^^ 
who obtained Tc ~ W/l^t, where W — 12t is the band- 
width for the 3-D Hubbard model, i.e. Tc/t sa 0.666. 



VIII. CONCLUSIONS 

In conclusion, we have investigated the ground state 
properties of the two-dimensional half-filled one band 
Hubbard model in the strong large-C/ to intermediate 
coupling limit i.e., away from the strict Heisenberg limit 
and antiferromagnetic phase diagram of the three dimen- 
sional Hubbard model using SU(2)xU(l) rotating ref- 
erence frame description. Our focus on systems in the 
strong to intermediate coupling regime, was motivated by 
the fact that weaker interactions are leading to increased 
electron mobility, which in turn should reduce the sta- 
bility of magnetic phases. Calculations with the Hamil- 
tonian for interacting electrons were reduced to calcula- 
tion of functional integrals with a phase-angular action. 
Collective bosonic fields are introduced by means of a 
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u/t 



FIG. 6: (Color online) Details of the phase diagram from Fig. 
[S] in the weak-coupling regime. Depicted is also the interac- 
tion dependence of the kinetic parameter v that enters the 
effective, J-induced hopping parameter tj. Note the vanish- 
ing of the true AF order for U/t ~ 0.676 and persistence of 
the gap Ac for arbitrary value of U /t. 



Hubbard- Stratonovich decoupling of the Hubbard inter- 
action and subsequent gauge transformation. Our imple- 
mentation for the Hubbard model is consistent with the 
spin rotation symmetry and simultaneously is able to re- 
produce the Hartrce-Fock result. One important techni- 
cal aspect arising in the construction of effective theories 
is that electron-defined operators in the bare high-energy 
theory are transformed into the composite particles sub- 
sequently employed in calculations within the effective 
low-energy theory. The inclusion of the quantum and 
spatial fluctuations has been shown to have a dramatic 
effect on transition temperatures and phase diagram. We 
have also compared the outcome of our calculations to a 
number of methods that were employed by other authors. 
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